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Abstract 

A "continuous measurement" Quantum Zeno Effect (QZE) in the context of trapped 
ions is predicted. We describe the physical system and study its exact time evolution 
showing the appearance of Zeno Phenomena. New indicators for the occurrence of 
QZE in oscillatory systems are proposed and carefully discussed. 
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Over the last few years the quantum dynamics of many models describing trapped 
ions have been investigated. The interest on this subject follows from the fact that 
such physical systems permit the observation of noteworthy quantum phenomena like, 
for example, quantum interference in Schrodinger cat states[l], decoherence[2,3,4] and 
quantum chaos[5,6]. In this letter our attention is focused on the quantum Zeno effect 
(QZE). In its standard form it consists of an inhibition of the system dynamics due to 
a continuous measurement process[7]. Misra and Sudarshan have in fact proved that 
in the limit of an infinite number of measurements performed in a fixed time interval, 
the system is frozen in its initial state[8]. Sometimes it happens that an appropriately 
large but finite number of measurements makes the initial state abandonment easier. 
In these cases we speak about Inverse Zeno Effect (IZE)[9,10]. 

Quite recently Facchi and Pascazio [9] have introduced a broader definition of Zeno 
effect which in this paper we refer to as Generalized QZE (GQZE). They consider the 
most general case in which an additional interaction, that may even not be a measure- 
ment interaction, is responsible for hindering the system dynamics. In order to better 
understand the physical meaning of GQZE, let us recall its definition more in detail. 
Suppose a physical system under scrutiny describable in terms of a hamiltonian model 



where K is a parameter and H meaB {K) is an interaction term which may even be as- 
sociated with a measurement process and such that Hm ea s(K = 0) = 0. By definition 
a GQZE occurs if, for K 7^ 0, there exists a time interval [tf , ] in which the survival 
probability is greater than the one correspondent to K = 0: 



H k = H + H, 



me as 
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P K (t) >P{t) 
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where 

P K (t) = | ty(0)| i>(t)) I 2 = I (^(0)| V(0)> I 2 (3) 

(Uk being the time evolution operator relative to H K ) is the probability to find the 
system in its initial state |\&(0)}, and 

P(t) = P*=°(t) (4) 

It is worth noting that this definition is applicable to both unstable and oscillatory 
systems. In these physical situations one can single out a Poincare time, T P (K), finite 
in the latter case and infinite in the former one. To claim the occurrence of a GQZE 
in oscillatory systems, in accordance with ref[9], one must require that — t%\ is of 
the same order of T P (K). 

In this paper we show the occurrence of a GQZE in the dynamics of a laser driven 
trapped ion. A formally equivalent result was found by Peres [11] in a different physical 
situation, while a standard QZE was brought to the light theoretically by Cook et al. 
and, experimentally, by Wineland et al.[12] in the context of trapped ions. 
Novel ways to characterize the appearance of Zeno Phenomena in oscillatory systems 
are moreover introduced and discussed in detail. 

Our physical system consists of a three-level ion confined in an isotropic three dimen- 
sional Paul trap interacting with a set of laser fields appropriately directed, polarized 
and tuned. 

In a Paul trap an inhomogeneous time dependent electric field and a quadrupolar 
static electric field confine a charged particle in a finite region of space wherein the 
ion moves subjected to an effectively three dimensional harmonic isotropic potential 
of frequency ujo [13,14]. 

The hamiltonian model describing such a situation is 

j—x,y.z i— 1,3 

where the annihilation (creation) operator aj(aj) describes the harmonic motion of the 
center of mass of the trapped ion along the axes j (j = x,y,z) and the i-ih internal 
level and its energy are denoted by \i) and hu>i respectively with i=l,2,3. 
A transition frequency u> between energy levels of Ho, involving the passage of the 
atom from the internal state \i) to the internal one \j) together with a simultaneously 
motional state modification is called m—ih red(blue) sideband when lu — (u>j — uii) = 
u) — ujij = mujo is negative(positive). If m = 0, that is when no bosonic transition is 
involved, to is called carrier frequency. 

Let us consider now the application on the trapped ion of a laser tuned to the first 
red sideband of the atomic transition 1-2 and another laser tuned to the first red 
sideband of the transition 2-3. The first laser is chosen of frequency co a = W12 — loq, 
wave vector k a , initial phase <f> a and directed along x and the second one, directed 
along y, characterized by uii = LO23 — wo, fcb, (j>b- The presence of the two lasers causes 
the appearance of a coupling between internal (electronic state) and external (center 
of mass) degrees of freedom of the trapped ion[15]. 

Supposing u) a ,Ub ^> wo, then the most important interaction term is the dipole one 
which in the Schrodinger picture can be cast in the form 

HI = -d-E= 

= ( E da |i>01)(£ E t e i «**-<* t+ *» + (h.c)) (6) 

i,j = l,3 Z — a,b 
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where E a and Eb describe amplitudes and polarizations of the applied fields, d is the 
atomic dipole operator and dij is the matrix element of d between |i) and 
In the interaction picture relative to the unperturbed hamiltonian Ho, expanding all 
the exponential operators and recalling that x(y) = e(a x ( y ) + we get 

j=0,oo 

after having discarded all the anti-rotating terms (Rotating Wave Approximation). 

In eq.(7) V a = ^ a ^ 12 , ^6 = r?i = fee with I = a, 6 and e = yQ^, M 

being the mass of the ion. The quantities ?? a , 776 measure the order of magnitude of the 
center of mass oscillation amplitudes with respect to the fields wavelengths. Finally 
ffij = |i) (j\ is the atomic transition operator. 

If Va,T)b «C 1 (Lamb-Dicke approximation), we are legitimated to truncate the series 
above up to the terms linear in r\ a and r\b, obtaining 

H{ = 7l o x |2> <1| + l2 a y |3> <2| + (h.c.) (8) 

where the initial phases <j) a and 4>b are chosen equal to \ to make real the quantities 

71 = — UQ, a e 2 rj a and 72 = —hQ, b e 2 7^. 

Implementing this hamiltonian model is certainly in the grasp of experimentalists. 
Generally speaking, one of the most attractive aspects of trapped ions systems is the 
relative easiness with which specific models describing the interaction between few 
level systems and few bosons may be realized, with respect to a lot of other physical 
situations like, for example, CQED[15]. 

In eq.(6) we have considered the sum of two laser fields and this leads to the effective 
hamiltonian model given by eq.(8). Increasing the number of laser fields and different 
choosing their frequencies, wavevectors, initial phases and amplitudes, we may real- 
ize effective hamiltonian models more and more complicated. For example, we may 
consider the physical situation wherein, in place of the laser 1-2 red sideband tuned 
previously used, we take two lasers tuned to the second 1-2 red sideband, it out of 
phase, with the same amplitudes and directed along the axis x' and y' rotated of J 
with respect to the axis x and y. The hamiltonian model describing such a physical 
system may be written down as H\ — ^\a x a y |2) (1| + yia y 3) (2| + (h.c.) [16]. 
On the basis of this example it is not difficult to convince oneself that it is always 
possible to engineer an appropriate laser beam configuration such that the vibronic 
coupling in the interaction picture may be represented in terms of the following general 
model 

Hi = 7l a " 1 2) <1| + 72 a ' |3> <2| + (h.c.) (9) 

where 

a V =a^a^a p z ' (P = r , 7) (10) 
and r , I are sets of three natural numbers. 

This is the hamiltonian model, on which we focus from now on. It generalizes eq.(8) 
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in the sense that more bosons and all the three vibrational modes are involved in the 
process. The related dynamical problem may be exactly solved for arbitrary r and I . 
Let us indicate the eigenstates of Ho with 



\n,aj = \n x ,n y ,n z ,cr) (11) 
where a = 1,2,3 refers to the atomic states. 

It is possible to see that in this basis the operator Hi is diagonal in blocks which are 
one dimensional if m — Ti < 0, two dimensional if rii — n > and rn — Vi — U < 0, 
and three dimensional if n 4 — n — h > with i = x,y, z. 

Here our attention is centered on a 3x3 representation of H\ having the form 



(12) 




in the subspace generated by | n , Ij , | n — r , 2j , | n — r — 1,3^ 
The non vanishing matrix elements appearing above are defined as follows 



a.(n,l|H 1 |n-?,2)= 7l / ft (13) 

y i = x,y,z 



{n-^,2\H A n-7-l,3)=,J ]J ^% (14) 

y i — x,y,z 



The parameters a and f3 measure the effective strength of the couplings between 
the three levels and are n-,r- and I -dependent. 

If Vi = k = for i = x,y, z, then the dynamical system described by the matrix (12), 
involves only the three electronic levels of the trapped ion. When, on the contrary, at 
least one of the laser is tuned to a sideband and the initial state of the center of mass 
is a Fock state such that n» — r» — h > for all i = x,y, z, then the hamiltonian matrix 
(12) describes the dynamics of a truly vibronic three-level system. 
It is straightforward to get the following time evolutions (in the interaction picture): 

\-> A , , , \B\' 2 + \a\' 2 coscut 1^ \ .a* . i_ _ \ 
|n, 1^> -> tp-, x (t) = it-J U I 71 ' 1 ) smujt | n - r ' 2 / + 

+ ^-(coswt-l)|n-r -T,3) (15) 

| n — 7, 2^> — > ip^_^ 2 (t) = coscjt | n — 7, 2^) — i% sinwt | n, 1^ + 

0* I _,_,-> \ 

— i— sinwt | n — r — I ,3j (16) 

|n - r - 7,3) ^_-j_ T]3 (t) = ^(coswt- 1) |n, l) + 

-j^ smut I n - r,2J + ^ | n - r - I ,3J (17) 
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where 



i=VW+W=\c\Vi + W (is) 



P = H TT ^ (20) 



with 

/? _ 72 TT »i - n! 

The quantity |x| expresses the relative strength of the coupling between the level 
In — r,2) with In — r — 1,3) and respectively and hence, as a and (3, is 



n-,r- and I -dependent and is determined both by laser parameters and by the initial 
condition. 

Such a property enables to obtain larger values of |x| with the same laser amplitudes 
only choosing different initial vibrational states, and hence allows to avoid possible 
difficulties related with the use of very intense laser fields. 

In passing we note that starting from the initial 3D zero point energy vibrational state, 
the Boulder group has succeeded in generating one-dimensional Fock states up to 16 
bosonic excitations. [15] 

In the limit |x| 3> 1, the time evolution of the state | n, 1^, as given by eq.(15), appears 
to be inhibited. This suggests us that when the coupling between levels 2 and 3 is 
greater than the one between 1 and 2, a GQZE occurs. 

In order to demonstrate the occurrence of such an effect in our system, in accordance 
with ref[9], let us compare our hamiltonian model given by eq.(9) with that given 

by eq.(l) making the correspondences H — > 71 a |2) (1| + (/i.e.), x * if ^ and 
H meas (K) -^K^a' 1 3) (2 1 + (h.c). 

Consider now the survival probability, P x (t), for different values of \, with |n, 1^ as 
initial state. Since the initial state is an eigenstate of Ho we have 

P*(t) = I (n, l| ^) | 2 = I (n, l| ^) | 2 (21) 

and therefore 

^)=(^+^) 2 (22) 

In view of eqs.(15)-(17) the time evolution of the trapped ion within the subspace 
under scrutiny, is governed by a single frequency w(x)- This naturally leads to define 
the Poincare-time, T p (x), of our system as the oscillation period, T(x), 

T P (X) = T(x) = = , = (23) 

Thus P x (0) = P x (T p (x)) = 1, which simply means that after a complete x-dependent 
Rabi cycle the system is found in its initial state. 

Figure 1 shows that for each value of |x > 1 there exists a time instant t x of the order 
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of T p (x) such that for te[0,t x ] we have P x (t) > P(t). Thus, at the light of ref[9], we 
may claim that our system manifests GQZE, under the condition previously discussed. 
It is worth noting that P x (t) exhibits an absolute minimum, m(x), in the time interval 
[0, T p (x)], which may easily be evaluated as 

fo if0<|x|<l, 



and occurs, for the first time, at 



arccos(-\ x \ 2 ) if0<|x|<l, 



•■W-T if M >i. (25) 

In figures 2 and 3 we plot the behaviours of m(x) and t m (x) respectively. 
The value of |x| at which m(x) begins to assume values different from zero is xl — 1 
and corresponds to the position of the maximum exhibited by t m {x)- It is interesting 
to note, in passing, that the non differentiability of t m (x) at Ixl — 1 originates from 
the coalescence of three stationary points (two minima and a maximum) of P x (t)[17]. 
We wish to emphasize that, as |x| grows up, the value of m(x) tends toward 1 and hence 
P x (t) tends toward 1 uniformly with respect to X- Such a property, namely m(x) — * 1 
as Ixl — » oo, describes in a very transparent way the effectiveness of the hindering 
process of the initial state time evolution as a function of x, and provides a new way 
to bring to the light the occurrence of a GQZE in physical systems characterized by a 
finite Poincare-time. 

It should be noted that this definition is stronger than the one given by Facchi and 
Pascazio[9]. This virtue proceeds from the simplicity of our physical system, which 
permits us to analyze in detail its dynamical properties and from the fact that we have 
restricted our GQZE-testing to oscillatory phenomena. 

Another way to construct an indicator of the GQZE-occurrence in finite Poincare-time 
systems is now proposed. Instead of considering the minimum of P x (t), we may use 
the mean value of the survival probability in a x-dependent Rabi cycle : 

r T(x i dt y 4 +'- 



This quantity has a very clear physical meaning representing the probability to find 
the system in its initial state independently of the time instant randomly chosen in 
[0,T( X )]. 

The behaviour of P(x) is plotted in fig. 4 and shows the existence of a minimum at 
x| = -^75 and the asymptotic tendency to 1 when |x| becomes larger and larger. This 
asymptotic behaviour may be reasonably taken as a clear indication that the system 
dynamics becomes more and more hindered "in average". In order to elucidate this 
statement, it is useful to introduce the quantity S(x, e) defined as the sum of ampli- 
tudes of the time intervals included in [0, T p (x)] wherein P x (t) < P(x) — £, e being a 
positive arbitrarily chosen x-independent real quantity. 

By "average hindering" we mean the tendency of the ratio to vanish as P(x) ~^ 1- 

Intuitively, since each P x {t) is upper-bounded by 1, if its temporal mean value, P{x), 

approaches this maximum, the function P x (t) tends to 1 almost anywhere. 

It is interesting to note that introducing, in analogy with eq.(26), the temporal mean 
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value of the probability P^ix) (-^Mx)) °f finding the trapped ion in the level n — r , 2 



( n — r — I ,3 J), the position of the absolute minimum of P(x) is just that in corre- 
spondence to which 



where the last equality follows from the normalization condition P(x)+P2{x)+Pa(x) = 
1. 

We wish to point out that the P(x)-test is weaker than the m(x)-test, in the sense 
that when m(x) grows up and approaches to 1, the parameter P(x) approaches to 1 
too but not vice versa. Nevertheless introducing P(x) appears to be useful in order 
to bring to the light a "global decrease" of P x (t) which is invisible when m(x)-test is 
considered. 

Summarizing, in this paper we have predicted that an ion confined in a Paul trap may 
exhibit a QZE when it is irradiated with appropriately configured laser beams. An 
important point that deserves to be emphasized is that our prediction concerns a con- 
tinuous measurements QZE which is different from the QZE "a la Misra-Sudarshan" 
predicted by Cook and experimentally found by Wineland[12]. The difference is that 
in our scheme only unitary evolutions are involved while in Cook's experiment de- 
scription projection operators are used and also the interaction with a series of short 
separated pulses is required. But the shorter is the distance between two pulses the 
more unfeasible is the series. On the contrary, a continuous interaction is simpler to 
realize and hence our scheme might be relatively easier to implement. 
The simplicity of the dynamical problem considered in this paper has given us the 
possibility of proposing novel indicators for the occurrence of the QZE. We feel that 
such novel ways of analysing the presence of QZE might be useful for different physical 
contexts provided that their dynamics is characterized by a finite recurrence time. 
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Figure 1: Survival probability P x (t) vs ^y-scaled time for different values of 
Ixl- Ixl = (black line), |x| = 1 (grey line), \x\ = 5 (grey dots), |x| = 10 (black 
dots) 
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Figure 2: Minimum m(%) of the survival probability vs |%| 
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